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Abstract 

We show that first approximations to the bulk viscosity 77,, are expressible in terms of factors 
that depend on the sound speed v^, the enthalpy, and the interaction (elastic and inelastic) cross 
section. The explicit dependence of 77,. on the factor Q ~ ^s) is demonstrated in the Chapman- 
Enskog approximation as well as the variational and relaxation time approaches. The interesting 
feature of bulk viscosity is that the dominant contributions at a given temperature arise from 
particles which are neither extremely nonrelativistic nor extremely relativistic. Numerical results 
for a model binary mixture are reported. 



1. Introduction 

Recent interest in the bulk viscosity rjy of strongly interacting matter stems from the ob- 
servation that in the phase transition from hadrons to strongly interacting quarks and gluons 77,, 
exhibits a drastic change f)^. The purpose of this work is to establish the dependence of t/i, on the 
sound speed v^, the thermodynamic properties of the system (particularly, the enthalpy) and the 
interaction (elastic and inelastic) cross sections between the constituents of a hadronic system. 

2. Bulk viscosity and the speed of sound 

2.1. Chapman- Enskog approximation ( Single component gas) 

In this approach, the first approximation to bulk viscosity can be written as 

= kT {a:l/2wf>) , (1) 

where 

«2 = -|z/.(r-3j + r|, z^—, h^— and r=-. (2) 

Above, h - h/c^ (given in terms of modified Bessel functions for a Boltzmann gas) is the reduced 
enthalpy and y is the ratio of specific heats at constant pressure and volume, respectively. The 
quantity wj, is the so-called omega integral which contains information about the cross section 
of the scattering particles. Explicitly (and for illustration, for elastic scattering). 



wf ^ —— \ di/rsinhVcosh'>/:/2zcosh(/r) t/0sin0 o-(i/',©) (1 - cos' 0) , (3) 
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where a-{ip, 0) is the differential cross section and j - | + | (-1)'; the others symbols are 



g = \{pi-P2) and P^{-p,p"fl^ (4) 
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sinhi/' = — and coshi^ = . (5) 
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using the relation 

dP 

r = 1 + ^ = 1 + V? . (6) 

Then, a2 in Eq. ([T) can be rewritten in terms of the speed of sound as 
Thus, the bulk viscosity takes the form 



ri.^kT — , (8) 



where a = -| {zh + l) and b - {zh - 4^. 



Limiting Situations 

It is instructive to consider the limiting cases of ultrarelativistic and nonrelativistic situations. 
For nearly massless particles or very high temperatures T, z - m/T «: 1. In this case, zh — > 4 
as z — > 0. Consequently, a — > - y and b — > 0. Utilizing these values, the bulk viscosity can 
be written as 

225 kT (I ,f 
'7v^^-^k-^? • (9) 



8 2w[,2H3 

Notice that for weakly interacting massless particles Vj |, so that 77,, — > 0. Note also the 
quadratic dependence of 77,, on v^. 

On the other hand, for masses such that z - m/T » 1, the coefficient a -jZ and b ^ -jz. 
These Umiting forms render the bulk viscosity as 



9|i-v?)\6(i-v?|.l 



(10) 



Weakly interacting, nonrelativistic particles are characterized by v, — » | so that, again 77,, — > 0. 

From the above analysis, we learn the important lesson that for a given temperature and 
for weakly interacting particles, intermediate mass particles contribute significantly to the bulk 
viscosity. It would be interesting to investigate the extent to which this conclusion is modified 
by strong interactions. 
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2.2. Variational and relaxation time approximations 

Here one starts from the general definition of the stress enegy tensor |!^ 



and dT - 



d^p 
(2;r)V" 



(11) 



where is the equilibrium distribution function, Vp is the velocity of a particle with momentum 
p and energy ep - In terms of the fluid velocity field m, the dissipative part of the energy 
stress tensor can be written as 



(12) 



Inserting the deviation function Sfp into the second part of the stress energy tensor and comparing 
Eqs. ( fTTT ) and ( fT2b . one gets |4] 



3/7 



2 2P 
p -m — 
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(13) 



where t is a momentum-independent relaxation time. In writing the rightmost equality, the 
identity v, - h/(cvT) and other integral identities have been used. Equation (T3[ lends itself to 
straightforward manipulations in the two limiting situations studied in the previous section. For 
m — > and — » |, one obtains the result that 77,, — > 0. In the case of massive particles, for which 
ep X m and — > |, the result 77,, ^ ^^^jA ^ is obtained. It is gratifying that the variational 
approach yields results similar to those obtained using the Chapman-Enskog approximation. 

It is easy to verify that similar conclusions emerge in the case that the expression for bulk 
viscosity features an energy dependent relaxation time Ta{ea) [1]: 



1 V r ^a(gfl) r/i 1 2\ 2 ly 



f 



(14) 



where the subscript a denotes the particle species. In conclusion, we learn that intermediate mass 
particles contribute the most to the bulk viscosity. 



3. Bulk viscosity of a binary mixture 

In the Chapman-Enskog approach, the first order approximation to the bulk viscosity of a 
binary mixture has the same formal expression as that for a one componet gas i3„JJ. Explicitly, 

[77,,]i=fer— , a2^xx— — \, (15) 
022 7(1) - 1 

where xi - ni/n, with n, denoting the number densities and y(,) denoting the ratio of specific 
heats of species /. The coefficient 022 differs from that in a one-component gas as it receives 
contributions from two sources: a22 = Ojj + Ojj . The quantity fl22 accounts for collisions 
between the same type of particles, and is calculated from the expressions given earlier for a 
one-component gas. The quantity accounts for collisions between different types of particles. 
Explicitly, 

16piP2 (1) , . 
«22 = >^1200(^12) ' (16) 
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where M - nii + 1112, pi - niini and the appropriate omega integral is 



W*2oo('^12) 



where P^^^ - rn^c^ + 
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c/0sin0i2Cri2(i/'i2,0l2)(l - COS ©12) ; 



rn^c^ + 2mim2C^ cosh 1^12 and ^12 = wic sinh 1^12 - m2C sinh 1/^12 



(17) 



Figure [T| shows illustrative results of bulk viscosity versus temperature for different mass 
configurations with constant cross sections. From the results shown, we learn that intermediate 
mass configurations contribute the most to the bulk viscosity. 



m, = 140 MeV 




T (MeV) 



Figure 1: Bulk viscosity versus temperature T in a binary mixture of dissimilar mass particles. The mass mi of particle 
1 is fixed at 140 MeV, whereas the mass m2 of particle 2 is varied as indicated alongside the different curves. Results are 
for a constant energy independent cross section of 1 fm^ for all collisions. 

Work is in progress to calculate the bulk viscosity of a hadronic mixture comprising of many 
hadronic resonances whose masses extend up to 2 GeV. 



Acknowledgments 

This research was supported by the Department of Energy under the grant DE-FG02-93ER40756. 



References 

[1] J. I. Kapusta, ArXiv :0809.3746v2. 

[2] W.A. Van Leeuwen et al, Physica A (1972) 65. 

[3] M. Prakash, M. Prakash, R. Venugopalan and G. Welke, Phys. Kept. 227, (1993) 321. 

[4] S. Gavin, Nucl. Phy. A 435, 826 (1985). 

[5] W.A. Van Leeuwen et al, Physica A XIV (1972) 233. 



4 



